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TECHNICAL MEMORANDUM 


the application of diffusion theory to the analysis of 

HYDROGEN DESORPTION DATA AT 25°C 


INTRODUCTION 


This study was directed toward the development of theoretical methods and their 
application to the explanation of experimental results concerning the uptake and 
elimination of gaseous hydrogen at ambient temperatures by high strength alloys. It 
is part of a broad program aimed toward an understanding of the nature of hydrogen 
embrittlement. Only the concentrations and behavior of the mobile or diffusable 
hydrogen were studied in the present work, the experimental information being 
obtained by electrochemical methods which are very sensitive to the small hydrogen 
concentrations obtained through charging at 25°C. Further work will include study 
°1 both -..-o mobile and trapped hydrogen, but, obviously, trapped hydrogen cannot 
tration* Sample depths S reater than those permitted by the overall hydrogen pene- 


Solutions of the Fick Law for diffusion in the X direction were employed for 
fw l™ pl * S AT With a11 oalcula tions being performed with computer programs written in 
the FOR IRAN language. Both the exact analytical solutions to the diffusion equa- 
tions, obtained as given by Barrer [1], and numerical methods, as used by the com- 
puter program PDEONE [2], were investigated. The analytical methods were 
employed because of their greater speed of computation and greater versatility 
Comparison of the analytical and numerical methods gave almost exactly the same 
numerical results. 


THEORETICAL METHODS 
Hydrogen Absorption 

All of the diffusion equations employed in this work have been described by 
Barrer m. When a slab is initially free of solute, the concentration at the faces are 
(at X - 0) and C 2 (at X = L), one has the example of diffusion of gas into a 

membrane of solid free of gas. For the present case, Cj = C 2 and the expression 
the series^* 1 hydrogen concentra tion in the sample after a charging time t is given by 
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where L is the sample thickness and D is the hydrogen diffusion coefficient. A plot 
of C versus t, with L = 0.232 cm and DM.3x 10~ 8 cm 2 /sec, with Cj = 1.0 ppm, 


is shown in Figure 1. Twenty-five iterations wore found to be sufficient for solu- 
tion of all diffusion equations. For concentration profiles, C(x) versus x, the follow- 
ing equation was used: 
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where S is the solubility of hydrogen in the metal at a given temperature and pres- 
sure and x is depth in the sample. The mean concentration is given by: 

L 

J' C(x)dx 


The concentration of hydrogen at the sample surfaces is immediately established at the 
solubility limit S for a given temperature and pressure or charging rate in the case 
of electrolytically charged samples. In practice the solubility S is obtained by divid- 
ing the mean (observed) hydrogen concentration by the value calculated by equation 
(1) for a given charge time with surfaces at unit concentration. The hydrogen con- 
centrations versus sample depth can then be obtained using equation (2), with the 
mean hydrogen concentration, as obtained using equation (3), being identically equal 
to the observed concentration. A plot of C(x) versus X is shown in Figure 2 for a 
Waspaloy sample charged for 1 hr at 5,000 psi (25°C). As seen from Figure 2, pene- 
tration into the sample is appreciable, the depth of penetration being almost 0.3 mm. 
The solubility obtained was 1.82 ppm, with the mean concentration, as obtained by 
equation (3) being 0.156 ppm, the observed hydrogen concentration. Thus, the type 
of curve predicted by diffusion theory is that illustrated in Figure 2, and will be 
referred to as a non-uniform distribution. A uniform distribution type, where the 
hydrogen concentration versus depth is constant, is illustrated in Figure 3, with a 
constant hydrogen concentration of 0.156 ppm. Such distributions would not be 
predicted by diffusion theory except after a very long charging time. 


Hydrogen Desorption 

Hydrogen desorption is appreciable at room temperature for both the nickel 
base alloys and the 4340 steel included in this study. It. is a quantity easily measured 
by sensitive methods, such as the EG&G-PARC corrosion measurement console, and 
provides important information concerning the hydrogen distribution in the metal 
after charging. For desorption problems, it is assumed that the concentration of 
hydrogen at the sample surfaces is zero for all time. The general equation for hydro- 
gen desorption is: 


cu,t) = * £ 
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where t is the desorption time in seconds. For a non-uniform initial distribution, or 
that predicted by diffusion theory, the vnlues of f(x') are those calculated by equa- 
tion (1) for a given charge time. The amount of hydrogen desorbed after a time t in 
given by: 


L 

Q(t) = Q h ^ - 0.0957341d J C(x,t) dx . (5) 

o 


Here, d is the metal density, and the conversion factor is that necessary to convert 

9 * 

from ppm-cm to coulombs /cm . Q HM is the total amount of hydrogen initially con- 
tained in the sample, and is measured experimentally. It is obtained as described by 
Zakroczymski [ 3] . For a uniform initial distribution , where the membrane contains a 
uniform concentration of hydrogen (C Q ) throughout, C(x,t) is given by: 
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and the total amount of hydrogen desorbed after time t is again given by equation 
(5). In practice, Q(t) versus t is measured experimentally through measurement of 
the current due to hydrogen desorption as a function of time, with Q(t) being 

obtained by integration of the current- time curve. Q H “ is the total amount of 

hydrogen desorbed after correction for cut-off at a finite time. The nature of the 
hydrogen distribution in the sample can be determined by fitting theoretical curves 
for non-uniform or uniform distributions, or a combination thereof, to the experi- 
mental curve. A typical case is illustrated in Figure 4 for 4340 steel charged elec- 

trolytically at 60 m.a/cm 2 for 1 hr. The curves are those calculated after a 10,000 sec 
desorption time. Curves are shown for non-uniform and uniform initial distribu- 
tions, with a hybrid curve representing a 59.2 percent uniform distribution, the 
observed quantity, also shown. 


EXPERIMENTAL 


The EG&G-PARC Model 350A corrosion measurement console was employed for 
all measurements in this work. The sample holder employed is shown in Figure 5. 

It presents a sample area of 1.0 cm 2 in contact with the solution in the electrolytic 
cell, in this case 0.1N NaOH. Samples of Waspaloy from two different sources were 
used. One sample was cut from a 1.27 cm (0.5 in.) diameter Waspaloy bolt on hand 
in the laboratory, while the other was cut from a used turbine disk, part S/N 2279525, 
obtained from a Space Shuttle Main Engine. The samples were essentially of the same 
composition and differed only in the method of preparation (heat treatment, etc.). 
Samples of Rene* 41 were used as obtained (annealed) without heat treatment. The 
4340 steel samples were heat treated to an ultimate tensile strength of 260,000 to 
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280,000 psi. All samples except that from the Waspnloy bolt were 1.59 cm (5/8 in.) 
in diameter, with nn average thickness of about 0.200 cm. 

The sample 3 were prepared by wet sanding with 150C and 220A grit silicon 
carbide paper and de-greasing in boiling trichloroethylene. Sample blanks were run 
at a constant potential of +0.25 V (NHE) in a 0.1N NuOH solution at 25°c. The period 
of measurement was 150,000 sec for each sample, with the value of the current being 
recorded at 500 sec intervals. 

Samples of all three materials were electrolytically charged with hydrogen in 

0.1N HgSO^ at two different charge rates, 60 ma/cm 2 and 1 ma/cm^. In addition, 

two samples of Waspaloy, one from the Waspaloy bolt described previously and the 
other from turbine disk part S/N 2279525, were charged with hydrogen gas at a 
pressure of 5,000 psi at 25°C. Data for the hydrogen-containing samples were 
collected in the same manner as those for the blanks, with the currents due to hydro- 
gen being obtained by subtraction. 

After each run, the current versus time data were read to a computer and 
saved. After data from both the hydrogen containing sample and its corresponding 
blank were obtained, the experimental curve, Q(t) versus time, was obtained by 
integration of the current- time curves corresponding to the current due to hydrogen 

only. Values of Q H “ , Q H ^, the observed concentration (C Q ) and the diffusion 

coefficient for hydrogen (D) were obtained according to the method developed by 
Zakroczymski [31, and described in a previous report [4]. 


RESULTS AND DISCUSSION 


The solubilities of hydrogen in the Rene' 41, Waspaloy and 4340 steel samples 
for electrolytically charged samples are shown in Table 1. As seen from the table, 
the solubility depends on the charge rate, being lower for the slow charge rates. 

This is analogous to the case for samples charged with high pressure hydrogen, 
where the solubility depends on the hydrogen pressure. The solubilities are about 
the same for the nickel-base alloys, Rene' 41 and Waspaloy, but are higher for a 
factor of about 4 for 4340 steel, which is known to be highly susceptible to hydrogen 
embrittlement [ 5] . The hydrogen distributions in the samples , both for high pressure 
charging and electrolytic charging, will be discussed separately for each metal. In 
all cases, theoretical curves for coulombs of hydrogen desorbed versus time will be 
shown for both uniform and non-uniform initial hydrogen distributions along with the 
curves for a hybrid distribution, if necessary, and the experimentally observed curve. 


Electrolytically Charged Rene' 41 

The desorption curves for a Rene' 41 sample charged for 4 hr at a rate of 

1 ma/cm 2 are shown in Figure 6 and those for a sample charged for 1 hr at 60 ma/cm 2 
in Figure 7. As seen from these figures, the hydrogen distribution in both cases is 
described very well by a non-uniform distribution which is in accord with diffusion 
theory . 


4 








Electrolytically Charged Wnspnloy 

1 he desorption curves for electrolytically charged Waspaloy at charge rates of 

t ma/cm and 60 ma/cm", respectively, are shown in Figures 8 and 9. As shown in 
Figure 8, fhe desorption curve for the slow charge rate is described very well by a 

theoretical curve corresponding to a 100 percent non-uniform distribution. For the 

o 

sample charged at 60 ma/cm , the experimental curve is best described by a hybrid 
curve, corresponding to a 10.2 percent uniform distribution, as shown in Figure 9. 
Thus, both experimental curves are essentially explained by diffusion theory, as was 
the case for the other nickel base alloy, Rene' 41. 


Waspaloy Charged 1 Hr at 5,000 psi 

The desorption curves for a Waspaloy sample cut from a used Space Shuttle 
Turbine Disk, part number S/N 2279525, are shown in Figure 10, while those for a 
sample cut from the Waspaloy bolt on hand in the laboratory are shown in Figure 11. 
The experimental curve for the turbine disk sample in Figure 10 is best described by 
a hybrid distribution which is 18.6 percent uniform, while that for the Waspaloy bolt 
sample in Figure 11 is completely explained by a 100 percent non-uniform distribution. 
Thus, the hydrogen distributions in Waspaloy, charged 1 hr at 5,000 psi at room 
temperature, are essentially in accord with those predicted by diffusion theory, as 
was the case for electrolytically charged samples. The solubility obtained for the 
Waspaloy bolt sample was 1.82 ppm, while that for the turbine disk sample was 1.28 
ppm. Thus, the hydrogen distributions and solubilities for the samples charged at 
5,000 psi are essentially the same as those obtained for electrolytically charged 

samples at a charge rate of 60 ma/cm 2 . 


Electrolytically Charged 4340 Steel. 

The hydrogen desorption curves for a 4340 steel sample charged 4 hr at 1 

ma/cm are shown in Figure 12. According to the results shown in Figure 12, the 
hydrogen distribution is 100 percent uniform in nature, a complete turnaround from 
the Rene' 41 and Waspaloy samples, and would not be predicted by the diffusion 
equations. While further study will be necessary to accurately explain this observa- 
tion, a possible explanation might be the following. 

For steady state flow, through unit area of any medium, the permeability is 
given by 


P = -D 


3C 
3 X 


(7) 


Thus, the permeability depends on the concentration gradient, as shown in equation 
(7), which would be quite high for 4340 steel where the solubilities are relatively 
large. Therefore, the absorbed hydrogen could pass rapidly through the surface 

layers into the interior of the sample, which would result in a more uniform distribu- 
tion. 
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The curves for a 4340 sample charged 1 hr at GO ma/cm are shown in Figure 
13. The observed curve is in agreement with a hybrid curve which is 59.2 percent 
uniform. Since a more uniform distribution would be predicted by diffusion theory 
for a longer charge time, the observations of Figures 12 and 13 are consistent in 
that respect. 


CONCLUSIONS 


The hydrogen desorption curves for Rene' 41 and Waspaloy samples, both 
nickel base alloys, are essentially described by non-uniform distributions and are 
in accord with those which would be predicted by diffusion theory for b'.th electro- 
ly tically charged samples and those charged at high pressure. The hydrogen solu- 
bility and distribution for a Waspaloy sample charged at 5,000 psi was essentially the 

a 

same as that obtained by electrolytic charging at 60 ma/cm . The hydrogen distribu- 
tion in 4340 steel charged electroly tically , was found to be uniform in nature, being 
100 percent uniform for a sample charged at a slow rate for a longer period of time. 
This situation would not be predicted from diffusion theory, and has been explained, 
for the time being, by the rapid diffusion of hydrogen across the surface layers, 
where the concentration gradient is very large in 4340 steel due to the much higher 
hydrogen solubility. This is analogous to the situation occurring in steady state flow. 
The hydrogen solubilities for all samples are higher, in the case of electroly tically 
charged samples, for high charging rates. This is analogous to high pressure charg- 
ing, where the solubility depends on the hydrogen pressure. 

From the results of this work, it is evident that the hydrogen desorption is 
completely explained by the nature of the hydrogen distributions in the samples. 

It was originally proposed [3] that the so-called "fast" hydrogen, or hydrogen which 
is desorbed at a faster rate than that which would be predicted by diffusion theory, 
is due to the formation of surface and subsurface hydrides which are unstable at 
room temperature. The hydrogen desorption rate was originally calculated assuming 
a uniform hydrogen distribution, which resulted in the discrepancy. Thus, the "fast" 
hydrogen is completely explained by the nature of the hydrogen distribution, and not 
by the formation of hydrides. 

Future work will include the study of cylindrical samples, with the determina- 
tion of both diffusable or mobile hydrogen and trapped hydrogen concentrations. 
Cylindrical samples will be used for a Leco Model RH-2 hydrogen analyzer and for 
the tensile tests to be performed on Satec Model C creep machines and, because of 
this, it will be necessary to study hydrogen desorption with the EG&G-PARC Model 
350A corrosion measurement console using cylindrical samples of the same diameter. 
Computer programs for cylindrical samples, which are analogous to those for flat 
samples, are now being developed. In addition, recourse might be made to hydrogen 
charging in molten salts, where very large hydrogen concentrations can be obtained 
at 150°C by charging in a eutectic mixture of NaHS0 4 .H 2 0 and KHSQ 4 [6]. 
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TABLE 1. OBSERVED SOLUBILITIES FOR ELECTROLYTICALLY 

CHARGED SAMPLES 
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Figure 12. Coulombs of desorbed versus time for eleetrolytically 
charged 4340 steel, slow charge. 
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charged 4340 steel, fast charge. 
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